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ABSTRACT 

The properties of the Galactic supernova remnant (SNR) SN 1006 are theoretically re-analyzed in 
the light of the recent H.E.S.S. results. Nonlinear kinetic theory is used to determine the momentum 
spectrum of cosmic rays (CRs) in space and time in the supernova remnant SN 1006. The physical pa- 
rameters of the model - proton injection rate, electron-to-proton ratio and downstream magnetic field 
strength - are determined through a fit of the result to the observed spatially-integrated synchrotron 
emission properties. The only remaining unknown astronomical parameter, the circumstellar gas num- 
ber density, is determined by a normalization of the amplitude of the gamma-ray flux to the observed 
amplitude. The bipolar morphology of both nonthermal X-ray and 7-ray emissions is explained by 
the preferential injection of suprathermal nuclei and subsequent magnetic field amplification in the 
quasi-parallel regions of the outer supernova shock. The above parameters provide an improved fit 
to all existing nonthermal emission data, including the TeV emission spectrum recently detected by 
H.E.S.S., with the circumstellar hydrogen gas number density Ah ~ 0.06 cm -3 close to values derived 
from observations of thermal X-rays. The hadronic and leptonic 7-ray emissions are of comparable 
strength. The overall energy of accelerated CRs at the present epoch is of the order of 5 % of the total 
hydrodynamic explosion energy, and is predicted to rise with time by a factor of ks 2. The relevance 
of CR escape from the SNR for the spectrum of the 7-ray emission is demonstrated. The sum of 
the results suggests that SN 1006 is a CR source with a high efficiency of nuclear CR production, as 
required for the Galactic CR sources, both in flux as well as in cutoff energy. 

Subject headings: acceleration of particles — cosmic rays — gamma-rays: general radiation 
mechanisms:non-thermal — shock waves — supernovae: individual(SN 1006) — 



1. INTRODUCTION 

The Galactic cosmic rays (CRs), with proton energies 
below a few 10 15 eV, are generally believed to be ac- 
celerated in shell- type supernova remnants (SNRs). To 
establish that the Galactic SNRs are indeed the main 
sources of the Galactic CRs one needs at least a hand- 
ful of SNRs with clearly determined astronomical pa- 
rameters, like the type of supernova explosion, the SNR 
age, the distance, and the properties of the circumstellar 
medium, together with the nonthermal emission charac- 
teristics from which the energetic particle spectra can 
be deduced with a theoretical model. The most suit- 
able objects are type la supernovae which have rather 
old and low-mass progenitors, not associated with star 
forming regions, and thus supposedly explode into a uni- 
form circumstellar medium - at least into one which is 
not modified by mass loss in the form of a stellar wind. 
The type of explosion also restricts the total explosion 
energy rather tightly and fixes the ejecta mass to a value 
close to the Chandrasekhar mass. 

One such object is the historical remnant SN 1006: 
the distance was determined us ing optical measurem ents 
with relatively high precision (jWinkler et al.l 120031 ) and 
all o ther astronomical paramete rs are quite well known 
(e.g. ICassam-Chenai et alj |2008). A special characteris- 
tic of SN 1006 is th e fact that it lies abou t 550 pc above 
the Galactic plane (Bocchi no et alJfeOllft in a very low- 
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density environment. This unusual property suggests 
that the ambient interstellar medium (ISM) is also free 
of density inhomg eneities from cooling instabilities (e.g. 
Ilnoue et al.ll2012ft . 

It has been shown earlier (jBerezhko et al.l I2009T) , in 
the following referred to as BKV09, that a nonlinear ki- 
netic theory of CR acceleratio n in SNRs (Berezhk o et al.1 
119961 iBerezhko fc VoHdll997ft 3 was consistent with all 
observational data available at the time. This included 
the H.E.S.S. energy flux density $(ei,e 2 ) = 2.5 x 
10 -12 erg cm -2 s" 1 of the Very High Energy (VHE:> 
100 GeV) 7-ray emission from one of the detected po- 
lar caps, integrated over the observed energy interval 
ei < e < e 2 , where e i = 0.2 TeV and e 2 = 40 TeV 
(N aumann-Godo et al.1 l2008f ) , as well as the fluxes of 
the nonthermal emission in the radio and X -ray regions 
measured with Cha ndra (| Allen et al.l 12004) and Suzaku 
(|Bamba et all 12008ft . Note that the above H.E.S.S. en- 
ergy flux density $(£1,62) was attributed in BKV09 to 
the whole remnant, whereas in fact it represents about 
half of the total flux. This is the reason why the val- 
ues of the ambient ISM density and of the SN explosion 
energy, estimated below, are correspondingly larger in 
comparison with BKV09. 

Compared with BKV09, in the following the physical 
parameters of SN 1006 - mainly the ambient gas density 
and the energy production efficiency - are re-estimated, 
based on the more recent measurements of the 7-ray 

3 Similar approaches were recently developed and applied to 
study the properties of other S NRs {ZTrakashvili fe Aharonianl 
[20101 ; FZirakashvili & Ptuskin 2011]). 
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spect rum and of the emission morphology (|Acero et al.l 
2010). In addition, the time profile of the acceleration 
efficieny is presented and the radial dependence of the 
7-ray brightness components is given. Finally the ra- 
dial profiles of the gas density and of the CR spectral 
energy density per unit logarithmic bandwidth are com- 
pared in a discussion of the energetic particle escape from 
the remnant. 

2. SUMMARY OF ASSUMPTIONS AND EARLIER RESULTS 

As a type la supernova SN 1006 presumably ejects 
roughly a Chandrasckhar mass M e j = 1.4M Q . Since the 
gas d ensity is observed to vary only mildly across the 
SNR (|Acero et al.1 12007). it appears reasonable to assume 
the circumstellar gas density and magnetic field to be 
roughly uniform. The ISM mass density po = m p NisM ~ 
l^mpA^H, is characterized by the hydrogen number den- 
sity ATh. As the most re liable distance estim ate the value 
d = 2.2 kpc is adopted (jWinkler et alJl2003T ). 

As in BKV09 a nonlinear kinetic theory of CR acceler- 
ation in SNRs (jBerezhko et al.1 1 19961 iBerezhko fc Volkl 
Il997| ) is applied. With one exception this theory in- 
cludes equations for the most important physical pro- 
cesses which influence CR acceleration and SNR dynam- 
ics: shock modification by CR backreaction, MHD wave 
damping and thus gas heating within the shock transi- 
tion, and synchrotron losses of CR electrons under the 
assumption of a strongly amplified magnetic field within 
the remnant B(t). 

The values of three scalar parameters in the govern- 
ing equations (proton injection rate 77, electron to pro- 
ton ratio K ep and the upstream magnetic field strength 
Bo) can be determined from a fit of the solutions that 
contain these parameters to the observed spatially inte- 
grated synchrotron emission data at the present epoch. 
Both, r\ and K op , are assumed to be independent of time. 
The parameter values for SN 1006, evaluated in this way, 
agree very well with the Chandra measurements of the X- 
ray synchrotron filaments and were obtaine d in BKV09 
by the analysis of the radio data compil ed by I Allen et al.l 
(|2004D (see also lAllen et al.ll200lL | 2008D and of the most 
accurate X-ray data of Cha ndra (| Allen et al.ll2004T ) and 
Suzaku (jBamba et alJl2008l ). 

3. RESULTS 

In order to explain the detailed 7-ray spectrum, val- 
ues of the hydrodynamic supernova explosion energy 
E sn = 2.4 x 10 51 erg and E sn = 1.9 x 10 51 erg are taken 
to fit the observed shock size i? s = 9.5 ± 0.35 pc and 
shock speed V s = 450 0±1300 km s" 1 (jMoffett et al.ll993t 
iKatsuda et al.l 120091 ) at the current epoch tsN ~ 10 3 yr 
for the ISM hydrogen number densities A^h = 0.08 cm -3 
and Nh = 0.05 cm -3 , respectively. These densities are 
consistent with the observed level of the VHE emission, 
as shown below. The best-fit value of the upstream 
magnetic field strength Bq = 30 //G is quite insensitive 
to TVh- The resulting current total shock compression 
ratios a for ATr = 0.08 cm -3 and A^h = 0.05 cm" 3 
become now a — 5.1 and 4.9, respectively, whereas 
the subshock compression ratios o~ s remain both close 
to a s = 3.7. Note, that the explosion energy value 
E sn = 2.4 x 10 51 erg is somewhat lower compared with 
E sn = 3 x 10 51 erg that one would expect extrapo- 
lating the value E sn = 3.8 x 10 51 erg obtained earlier 
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Fig. 1. — Overall energy E c of accelerated CRs, normalized 
to the total hydrodynamic energy release E Bn , as a function of 
time, calculated for Nn = 0.08 cm -3 (solid line) and for TVh = 
0.05 cm -3 (dashed line). Here and in all following figures the 
quantity Bq = B c \/<r denotes the amplified upstream magnetic 
field strength, where a is the overall shock compression ratio and 
_B<j is the downstream magnetic field strength. The vertical dotted 
line shows the current evolutionary stage tsN- 

(|Ksenofontov et all 120051 ) for A^h = 0.1 cm" 3 according 
to the expected dependence E sn ex Ah- However this 
difference is well within the range determined by the un- 
certenties of the observed values R s and V s taking into 

1/5 

account the relation R s oc £ s „ . 

Since the properties of the accelerated CR nuclear and 
electron spectra and their dependence on the relevant 
physical parameters, as well as the dynamical properties 
of the system were described in detail in BKV09, they 
will not be discussed here once more. 

3.1. Acceleration efficiency 

In Fig. [T] the time dependence of the fractional energy 
E c /E sn contained in accelerated CRs during the SNR 
evolution is presented. Note that the value of E c is re- 
duced by a factor of / rc = 0.2 (see below) compared 
with the value calculated within the spherically symmet- 
ric model. According to Fig. [TJ E c /E sn » 0.05 and 
0.065 for Af H = 0.05 cm" 3 and Af H = 0.08 cm" 3 , re- 
spectively. This is lower than the value « 0.1 which is 
required on average for each SNR for their population to 
be the main source of CRs in the Galaxy. The reason is 
that SN 1006 is a quite young object in an evolutionary 
sense: according to Fig. [T] E c (t)/E sn is expected to ap- 
proach this canonical value in the subsequent evolution. 
On the other hand, the requirement for an average effi- 
ciency of w 0.1 is based on an assumed Galactic average 
value E sn = 10 51 erg which is about one half of our value 
E sn = 2a:10 51 erg. Therefore the calculated efficiency 
E c /E sn w 0.05 fulfills the average requirement even for 
the present epoch. 

3.2. Overall nonthermal spectra 

Fig. [5] illustrates the consistency of the synchrotron 
and 7-ray spectra, calculated with the best set of pa- 
rameters ( 7] = 2 x 10" 4 as well as K cp = 4.5 x 10" 4 
for Nft = 0.05 cm" 3 , and r\ = 2 x 10" 4 as well as 
K cp = 3.2 x 10" 4 for A^ = 0.08 cm" 3 ), with the ob- 
served spatially in tegrated spectra. The H.E.S.S. data 
(|Acero et all l2010f ) for the NE and SW limbs, respec- 
tively, have been multiplied by a factor of 2, in order 
facilitate comparison with the full deduced 7-ray flux. 
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As can be seen from Fig. [2j the calculated synchrotron 
spectrum fits the observations both in the radio an d the 
X-ray ranges (|AUen et al.ll2004t Ba mba et al.ll2008D very 
well. Note that the Chandra flux (|Allen et all 12001) is 
from a small region of the bright NE limb with mini- 
mal contributions from thermal X-rays; t his X-ray flux 
was n ormalized to the overall Suzaku flux ([Bamba et al.l 
2008) at energies e > 2 keV (for details, see BKV09). 

The only important parameter which can not be de- 
termined from the analysis of the synchrotron emission 
data is the external gas number density Nn: Fig.[2]shows 
that the spectrum of synchrotron emission is almost non- 
sensitive to the ambient gas density. Consequently, nu- 
merical solutions have been calculated for the pair of val- 
ues JVh = 0.08 cm~ 3 and jVh = 0.05 cm -3 which appear 
to bracket the density range consistent with the H.E.S.S. 
7-ray measurements. It should be noted that this es- 
timate for N-r is by a factor of w 1.5 larger than the 
estimate of BKV09. The reason is that in the analysis of 
BKV09 the total calculated 7-ray energy flux was com- 
pared with the energy flux from the NE rim. However, 
the latter flux corresponded only to about 50% of the to- 
tal observed energy flu x, as it had been reported for the 
H.E.S.S. instrument by iNaumann-Godo et ail (|2008l ). A 
detailed comparison of the 7-ray spectrum with the cal- 
culated spectrum will be made in the context of Fig. 



3.3. Gamma-ray morphology 

The 7-ra y morphology as fo und in the H.E.S.S. mea- 
surements (|Acero et all 120 lOT) . is consistent with the 
prediction of a polar cap geometry on account of the 
strongly preferred injection of nucle ar particles in th e 
quasi-parallel regions of the shock (iVolk et al.l l2003h . 
Such a geometry has also been found experimentally from 
an analysis of the syn c hrotro n mor phology in hard X-rays 
by iRothenflug etail ([2004) and ICassam-Chena'i et al.l 
(200|). This means that the 7-ray emission calculated 
in the spherically symmetric model m ust be renormalized 
(reduc ed) by a factor / re rs 0.2, as in iKsenofontov et al.l 
(2005) and BKV09. This renormalization factor is ap- 
plied here. The total 7-ray flux from the source is the 
sum of the fluxes from these two regions. 

This morphology is also a key argument for the exis- 
tence of an energetically dominant nuclear CR compo- 
nent in SN 1006, because only such a component can 
amplify the magnetic field to the observed degree. The 
accelerated electrons are unable to amplify the magnetic 
field to the required level (BKV09). 

The question, whether these bright NE and SW re- 
gions of SN 1006 represent quasi-par allel portions of the 
SN shock or not, is still debated JlPetruk et al.l 120091: 
ISchneiter et al.ll2010t iMorlino et al.ll2010[ ). On the other 
hand iBocchino et al.l ( 20111) have recently argued from 
the radio morphology that the radio limbs - morpho- 
logically similar to the X-ray and 7-ray limbs - are po- 
lar caps and that electrons are accelerated with quasi- 
parallel injection efficieny Like the aforementioned 
theoretical arguments and X-ray analyses, this argues 
against a scenario, where t he magnetic field is perpendic - 
ular to the sho c k nor mal ([Fulbright fc Reynolds! Il990f ) . 
IBocchino et al.l (|2011f) also concluded that the remain- 
ing asymetries (converging limbs and different surface 
brightness), clearly visible also in the H.E.S.S. 7-ray data 



(|Acero et al.ll2010l ), could be explained by a gradient of 
the ambient ISM magnetic field strength. 

Additional arguments which confirm that the NE and 
SW regions correspond to quasi- parallel shock regions 
are given by iRevnoso et al.l (|201 lh . Their analysis of the 
polarization of radio emission led them to the conclusion 
that the magnetic field in SN 1006 is radial at the NE 
and SW lobes but tangential in the SE of the radio shell. 
In addition they established the maximum fractional po- 
larization in the SE which implies that the magnetic field 
is highly ordered there. On the other hand the low frac- 
tional polarization in the lobes suggests a considerable 
randomization of the magnetic field. The fact that a 
turbulent magnetic field coexists with the brightest syn- 
chrotron emission in the SW strongly supports efficient 
acceleration of the CR nuclear component, followed by 
magnetic field amplification which provides its own di- 
rectional randomization. 

3.4. Relative contributions to the gamma-ray flux 

Fig. [3] shows the 7r°-decay, the inverse Compton (IC), 
and the total (7r°-decay plus IC) 7-ray energy spectra 
of the remnant, calculated for TVjj = 0.05 cm~ 3 and 
Ah = 0.08 cm~ 3 , respectively. It can be seen that the 
7-ray spectrum produced by the nuclear CRs is rather 
close to the IC emission spectrum produced by CR elec- 
trons alone, especially in the case of Ah = 0.05 cm~ 3 . 
Even for the case Ah = 0.035 cm -3 which is close to 
the minimum density Ah = 0.03 cm -3 for the remnant, 
as argued by lAcero et al.l (|2007( ). the total predicted 7- 
ray flux is still quite similar in spectral shape to the 
total flux for Ah = 0.05 cm -3 (and essentially also for 
A H = 0.08 cm" 3 ), as shown in Fig.6 of BKV09. Un- 
less the actual gas density would be small compared to 
Ah = 0.03 cm -3 , it would therefore be very difficult 
to observationally distinguish a hypothetical dominant 
IC scenario from the mixed 7-ray spectrum discussed 
here. Therefore, in the VHE range - except perhaps in 
the last decade of energy, near the cutoff, see subsec- 
tion 3.6 - the observed 7-ray spectrum alone is not able 
to discriminate between the hadronic 7r°-decay and the 
leptonic I C 7-ray components. How ever, it was already 
shown by Ksenofon tov et aD (|2005h that such a low to- 
tal VHE emission flux, with a highly depressed IC 7-ray 
flux due to the synchrotron losses of high energy elec- 
trons, is only possible if the nuclear CR component is 
efficiently produced with accompanying strong magnetic 
field amplification. 

3.5. Radial profile of gamma-ray brightness 

A possibility for an experimental discrimination be- 
tween 7r°-decay and IC 7-rays is given by the measure- 
ment of the radial p rofile of the 7-ray b rightness. It 
was shown early on ([Berezhko et al.|[2002f) that the ra- 
dial profile of the TeV 7r°-decay 7-ray emission has a 
peak near the remnant rim with a width of about 20% of 
the remnant radius. On the other hand, the peak of the 
IC 7-ray emission was predicted to be much narrower, 
actually narrower by an order of magnitude. This lat- 
ter characteristic is physically identical to the filamen- 
tary structu re detected subsequently with Chan dra in 
keV X-rays dBamba et al.l 120031: lLong et al.l 120031 (and 
for Cas A bv lVink fc Laming] (|2003l )). because the keV 
synchrotron emission and the TeV 7-ray IC emission are 
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Fig. 2. — Spatially integrated spectral energy distribution of SN 1006, calculated for Nn = 0.08 cm -3 (so lid line) and iVjj = 0.05 cm -3 
(dashed line), and compared with observational data. The observed radio spectrum is from a compilation of Allen ct al. (2004). The form 
of the Chandra X-ray spectrum w as measured for a sm all region of the bright northeastern (NE) rim of SN 1006 < Allen et al.H2004T) , whereas 
the Suzaku spectrum (blue color) (Bamba ct al. 2008) is for the entire remnant. The corresponding Chandra X-ray flux has been mu l tiplied 
by a factor of 60 (red colour) such as to be consistent with the Suzaku flux for energies e > 2 keV. The H.E.S.S. datafAccro ct al. 1 120101 ) 
are for the NE region (green colour) and the SW region (magenta colour). Both these H.E.S.S. data sets have been multiplied by a factor 
of 2, in order to facilitate a comparison with the global theoretical spectra. 

agrees with the H.E.S.S. observations for p > 0.5i? s but 
exceeds the H.E.S.S. values considerably at low radial 
distances p < 0.5i? s . The reason for this discrepancy is 
simple. As it has been said in subsection 3.3, only those 
nuclear CRs that occupy two cones with opening angle 
20° , with symmetry axes going from the remnant center 
toward the NE and SW directions, respectively, should 
be taken into account to model the hadronic emission of 
the actual remnant due to the strong angular dependence 
of the nuclear injection rate. The radial profiles of the 
7-ray emission lobes, corresponding to the NE-SW di- 
rections produced by these CRs, are well consistent with 
the H.E.S.S. data (Fig. [4J . In this case the line of sight at 
p < 0.7R S does not intersect the outer spherical region of 
the remnant which contains most of the accelerated CRs 
that provide the brightness depression at p < 0.7 R s . 

A discrimination between these 7r°-decay and IC 7-ray 
components would require a 7-ray instrument with an 
angular resolution that is an order of magnitude higher 
than that of H.E.S.S., or of any other existing 7-ray 
instrument. It is, however to be noted that the measured 
radial H.E.S.S. profile, which gives a width of the mea- 
sured peaks of 20% of the remnant radius (jAcero et al.l 
2010), is clear evidence that the nuclear CR component 
is accelerated efficiently. Indeed, if the energy content 
in CR nuclei were small due to inefficient proton accel- 
eration one would have to conclude that the magnetic 
field in the remnant is not amplified to any substan- 
tial degree. In such an inefficient scenario electrons with 
energies of tens of TeV, which produce the TeV 7-ray 
IC emission, would be extremely smoothly distributed 




Fig. 3. — Total (7r°-decay + IC) (black lines), 7r°-decay (red lines) 
and IC (green lines) differential 7-ray energy fluxes as a function 
of 7-ray energy, calculated for the ISM hydrogen number densities 
Nn = 0.05 cm~ 3 (dashed lines) and jVjj = 0.08 cm~ 3 (solid lines), 
respectively, and for the injection parameters derived from the fit 
of the overall synchrotron spectrum. 

produced by the same electrons, with energies approach- 
ing 100 TeV. 

In Fig. |4] radial profiles of the 7-ray brightness 
Jj(p) are shown as functions of radial projection distance 
p corresponding to energies e 7 > 0.5 TeV, calculated for 
Ah = 0.05 cm~ 3 . It is seen that the profile of the IC com- 
ponent is indeed dramatically thinner than the profile of 
the 7r°-decay component. The profile of the total 7-ray 
emission, which is the sum of these two components, has 
a width at half maximum of about 0.15i? s that is smaller 
than that observed by H.E.S.S. The total radial emission 
profile, smoothed to the H.E.S.S. point spread function, 



across the remn ant - almost uniformly in th e spherically 
symmetric case (|Berezhko et al J 12002, 12001 . Therefore 
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Fig. 4. — Radial dependence of the7-ray brightness for 7-ray 
energies e 7 > 0.5 TeV. IC (green color) and 7r°-decay (red color) 7- 
ray components together with their sum are shown by thin dotted, 
thin dash-dotted and thin dashed lines respectively. These radial 
profiles, smoothed to the H.E.S.S. point spread function, are shown 
by corresponding thick lines. The thick solid line corresponds to 
the total 7-ray emission that originates i n two cone regions . The 
H.E.S.S. measurements are shown as well (Accro ct al. 2010). The 
angular radius 0.24° corresponding to the position of the peak of 
the H.E.S.S. 7-ray profile is taken equal to p = 0.85R S - 

the expected radial width of the 7-ray emission region 
would be considerably larger than that expected from 
efficiently accelerated protons. 

3.6. External density 

In order to exhibit the most detailed comparison of the 
calculated 7-ray spectra with the observations, in Fig. [5] 
the theoretical spectra for two ambient densities are pre- 
sented together with the observed H.E.S.S. spectra for 
each of the polar caps on an expanded energy scale. For 
ease of comparison, the theoretical spectra are divided 
by a factor of two, which would correspond to the emis- 
sion from a single polar cap region for an ideal dipolar 
morphology. 

The H.E.S.S. spectra from the two regions are compat- 
ible with decreasing power-law distributions dF^/de-y oc 
e~ r , with T ~ 2.3 and somewhat different fluxes <&(> 
ITeV), in the detected energy range 11.4 < loge^< 13.2 
(|Acero et al.ll2010t ). As Fig. El Fig. M and Fig. U show, 
the H.E.S.S. data are also compatible with the more com- 
plex total theoretical 7-ray spectral energy flux density 
over this range that has a maximum at loge 7 ~ 12.6. 
Above 10 GeV the theoretical spectral energy flux can 
be approximated analytically by a power law with an 
exponential cutoff: 



dFj/de^ oc e exp (- 



yAc), 



with cutoff energies e c = 30 TeV and 37 TeV for N H = 
0.05 cm -3 and Ah = 0.08 cm -3 , respectively. These two 
analytical curves represent the best fit to the H.E.S.S. 
data. At the same time they coincide within 4%, respec- 
tively, with the theoretical curves presented in Fig. O 
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Fig. 5. — Spatially integrated 7-ray energy flux from SN 1006, 
divided by a factor of 2, calculated for the two ISM hydrogen num- 
ber densities Nn = 0.05 cm -3 (dashed lines) and ./Vh = 0-08 cm " 3 
(solid lines), respectively. The H.E.S.S. data (Accro ct al. 2010) for 
the NE region (green colour) and the SW region (magenta colour) 
are shown as well. 

The quality of fit to the H.E.S.S. data is characterized 
by the values y^/dof = 1.05 and 2.2 for the NE and 
SW part of the remnant respectively. It is comparable 
with the quality of fit by the pure power law distribu- 
tions dF^/dej oc e~ r , with T ~ 2.3 and different fluxes 
(jAcero et al.1 |2010T l which have, according to our cal- 
culation, a x 2 /dof — 1.3 and 2.5 for the two regions, 
respectively. 

In order to discriminate between those different spec- 
tral forms the observational energy range should be ex- 
tended both below the present H.E.S.S. threshold en- 
ergy and above the present maximum energy detected 
by H.E.S.S. 

According to Fig. [5] the H.E.S.S. data are consistent 
with an ISM number density in the interval 0.05 < -/Vh < 

0. 08 cm~ 3 . There is also some indication that the ac- 
tual number density is closer to the lower end of this 
interval, say A^h ~ 0.06 cm~ 3 . However this is a weak 
conclusion, given the approximations made in the the- 
oretical model in terms of the reduction factor f re (see 
above). It should be noted nevertheless that the lower 
end of this interval is also preferred from the point of 
view of supernova explosion theory: the corresponding 
hydrodynamic explosion energy E sn « 1.9 x 10 51 erg is 
close to the upper end of the typical range of type la 
SN explosion energies that vary by a factor of about two 
(jGamezo et al.ll2005t iBlinnikov etal1|2006D. In ad dition, 
from their X-ray measurements also lAcero et al.l (|2007l ) 
concluded that a value Ah ~ 0.05 cm -3 is probably rep- 
resentative for the ambient density around SN 1006. 

A solid observational conclusion is that the total emis- 
sion from the NE limb exceeds that from the SW limb 
(jAcero et al.ll20l6T ) . Given that the hadronic 7-ray flux is 
proportional to A^, the result in Fig. 5 suggests that this 
excess is due to a gas density difference SN-r ~ 0.02 cm~ 3 
between the two 7-ray emission regions across the SNR, 

1. e. a relative difference of ~ 30%. Given that the den- 
sity difference between the thermally emitting NW and 
SE regions is of the order of 0.1 cm~ 3 (jAcero et al.ll2007h 
such a minor density difference is indeed quite plausi- 
ble. However, it does not explain the deviations from 
full dipolar symmetry in terms of a noticeable conver- 
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0.9 1 1.1 

Fig. 6. — Radial profiles of the gas number density N s (in units 
of the postshock number density N g 2 ) and of the proton distribu- 
tion / in momentum p, multiplied by e 4 , for three values of the 
proton energy e, in eV. The figure represents the actual depen- 
dence of the distribution function on energy, whereas the common 
normalization is arbitrary. 

gence of the limbs and thus should simply complement 
the gradient in the m agnetic field strength claimed by 
iBocchino et al.l (|2011l ) from radio observations. 

3.7. Escape 

It should be noted that the cutoff of the theoretical 
7-ray spectrum (Fig. [5]) begins already at e 7 rs 10 13 eV. 
This value is unexpectedly low since the overall pro- 
ton power-law spectrum extends up to a proton energy 
e rs 10 15 eV (BKV09). The reason is the progressive 
decrease of overlap between the radial gas density pro- 
file iVg(r) = p(r)/m p and the radial profile of the pro- 
ton distribution function f(r, e) at energies e > 10 14 eV. 
As Fig. |6] shows, the gas number density has a peak 
Ng2 = ctNism a t the subshock position r = R s and de- 
creases towards the remnant center by a factor of 2 at a 
distance l g rs 0.06i? s from the subshock position. (The 
thickness of the thermal gas precursor is about 10" 4 i? s ; 
therefore it is hardly visible in Fig. [BJ) The correspond- 
ing radial dependence of the proton distribution in mo- 
mentum p, multiplied by e 4 , i.e. of the spectral energy 
density cp 4 f(r,p), with energies e rs cp < 10 14 eV, has 
a similar behaviour; it shows a relatively small particle 
number in the upstream region r > R s . At higher ener- 
gies e > 10 14 eV the thickness of the radial distribution 
increases with energy, so that for e rs 10 15 eV its overlap 
with the gas density profile is by a factor of two lower 
than for e = 10 14 eV. This leads to a lower cutoff of the 
7-ray spectrum compared with what one would expect 
in the simple case of energy-independent overlap of CR 
protons with the gas distribution. 

This spread of CRs into the upstream region r > i? s , 
which for the highest energies e ~ e max (i) is faster than 
the shock expansion Rs(t), represents the diffusive es- 
cape of CRs from the expanding SN R and was pre- 
dicted by iBerezhko fc Krvmskvl (|1988[ ). Such a CR es- 
cape is expected in particular in the Sedov phase when 
the maximum energy e max (i) oc BqR s V s of CRs, accel- 
erated in the given evolutionary phase, decreases with 
time, because t he value of e max is shock-size limited (see 
lBerezhkolll996i for details) rather than time limited. At 
any given phase t > to (where to is the sweep-up time) 
efficient CR acceleration at the SN shock takes place 
only for energies e < e max (i), whereas for CRs with en- 



ergies e max (i) < e < e max (to), produced during earlier 
times, the acceleration process becomes inefficient and 
these CRs expand into the upstream region. Within the 
SNR model discussed in this paper CR escape is rela- 
tively slow because the model assumes Bohm diffusion 
in the amplified magnetic field Bo everywhere upstream 
for any CR energy. In reality this is expected to be true 
only for CRs with energies e < e max (i) in the vicinity of 
the shock, where these CRs produce significant magnetic 
field amplification. At large distances upstream of the 
shock, r — R s ^> 0.1R S , or/and for higher CR energies 
e > e max (t), CR diffusion is much faster than Bohm diffu- 
sion. Therefore in actual SNRs CR escape is expected to 
be faster and more intense than the present model pre- 
dicts (see iPtuskin fc Zirakashyihl 120031: lDrurvll2011l for 
more detailed considerations). Since SN 1006 is only at 
the very beginning of the Sedov phase this underestimate 
of the magnitude of escape is however not very critical. 

4. SUMMARY 

Based on a nonlinear kinetic model for CR acceler- 
ation in SNRs the physical properties of the remnant 
SN 1006 were examined in a detail which corresponds 
to that of the available observations. Since the relevant 
astronomical parameters as well as the synchrotron spec- 
trum of SN 1006 are measured with high accuracy, the 
values of the relevant physical parameters of the model 
can be estimated with similar accuracy for this SNR: 
proton injection rate r\ rs 2 x 10~ 4 , electron-to-proton 
ratio K ep rs 3.8 x 10 -4 and downstream magnetic field 
strength rs 150 ^G. 

As a result the flux of TeV emission detected by 
H.E.S.S. is predicted to agree with an ISM hydrogen 
number density Ah ~ 0.06 cm~ 3 , consistent with the ex- 
pectation from X-ray measurements. The corresponding 
hydrodynamic SN explosion energy _E sn w 1.9 x 10 51 erg 
is somewhat above but fairly close to the upper end 
E sn = 1.6 x 10 51 erg of the typical range of type la 
SN explosion energies. The efficiency of CR production 
is predicted to be between 5 and 6 percent up to the 
present epoch and expected to approach of the order of 
10 precent during the further evolution. Normalized to 
a standard value E sn — 10 51 erg, already the present 
efficiency is about 10 percent. 

The magnetic field amplification properties of this SNR 
can be understood as the result of azimuthal variations of 
the nuclear ion injection rate over the projected SNR cir- 
cumference in terms of injection at quasi-parallel shocks 
only and the corresponding acceleration. This predicts 
the dipolar 7-ray emission morphology and is compati- 
ble with a polar cap- type X-ray synchrotron morphology. 
The 7-ray morphology is the key argument for the ex- 
istence of an energetically dominant nuclear CR compo- 
nent in SN 1006. The magnetic field amplification allows 
the proton spectrum to extend to energies ~ 10 15 eV. 
The difference between the 7-ray emission from the NE 
polar region and that from the SW polar region can be 
attributed to a small density difference between those 
regions. 

It is shown that the high-energy tail of the 7-ray 
spectrum is affected by CR escape from the SNR inte- 
rior, consistent with the comparatively low 7-ray cutoff 
near 10 14 eV. 

It is also concluded that the radial profile of the TeV 
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7-ray emission measured by H.E.S.S. is evidence that the 
nuclear CR component is indeed efficiently produced. In 
the opposite case of inefficient nuclear CR production 
the magnetic field would not be expected to be amplified 
and the radial profile of the IC-dominated 7-ray emission 
would be expected to be significantly smoother than ob- 
served. 

The sum of all these results suggests the conclusion 
that SN 1006 is a CR source with a high efficiency of 
nuclear CR production, required for the Galactic CR 
sources, both in flux as well as in cutoff energy. 
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